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Abstract 

Using the self-consistent Born approximation, and the corresponding wave 
function of the magnetic polaron, we calculate the quasiparticle weight corre- 



a 
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sponding to destruction of a real electron (in contrast to creation of a spinless 

holon) , as a funtion of wave vector for one hole in a generalized t — J model 

> 

<0 and the strong coupling limit of a generalized Hubbard model. The results 

O ' 

are in excellent agreement with those obtained by exact diagonalization of a 

o 

£-~. ' sufficiently large cluster. Only the Hubbard weigth compares very well with 

photoemission measurements in S^CuC^C^. 
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The problem of a single hole in an antiferromagnetic background has been a subject of 
considerable interest since the discovery of high-T c systems . One of the most powerful tools 
for this study is the self-consistent Born approximation (SCBA) [|1]-||]. Excellent agreement 
has been obtained between the position of the lowest pole of the holon Green function of the 
SCBA and the quasiparticle dispersion obtained by exact diagonalization of small systems 
H|-[5|]. An important advance in the understanding of the SCBA has been the explicit 
construction of the corresponding wave function by Reiter ||. 

The interest on the problem has been revived by recent angle-resolved photoemission 
experiments on insulating Sr2Cu02Cl2, in which the hole dispersion and quasiparticle weight 
have been measured 0. While it was clear that the "bare" t — J model was unable to explain 
the observed dispersion, several works have appeared fitting the experimental dispersion 
using generalized t — J models |5|,@,[|, a generalized Hubbard model |T(J and the spin- 



fermion (or Kondo-Heisenberg) model for the cuprates JTTJ. Except for the fact that the 



band width is ~ 10% narrower than the experimental result if the experimental value of 
J is taken ]T2|, the generalized t—J model including hopping to second and third NN 



and the three-site term t", reproduces well the experimental dispersion [|5],|9| and also other 
properties of the spin-fermion and three-band Hubbard models |]T2[. A consistent picture 
of the observed spin and charge excitations has been obtained using a generalized one-band 



Hubbard model 10 



However, very little attention has been devoted to the explanation of the intensity of the 
observed quasiparticle peaks. This task is difficult for the following reasons: i) exact results 
for quasiparticle intensities in sufficiently large clusters (containing more than 16 unit cells, 
as discussed below) exist only for the "bare "t — J model and only at a few wave vectors, 
ii) The SCBA provides the Green function of the spinless holon, while the Green function 
of the real particles contain spin-wave excitations and simple decoupling approximations do 
not provide reasonable results. The holon weights are the same for wave vectors differing in 
(it, 7i") contrary to experiment, iii) While a lot of work has been devoted to the mapping of the 
three-band Hubbard model for the cuprates to low-energy effective models, less attention has 



been devoted to the mapping of the corresponding operators p2| , ^4[p!5[ . This information 
as well as the photoionization cross sections for Cu and O are necessary if accurate weights 
are wished. 

In this paper we calculate the photoemission quasiparticle weight for removing an elec- 
tron, as a function of wave vector in generalized t — J and strong-coupling Hubbard models, 
using the SCBA and the wave function of the polaron ||. The Hamiltonian has the form 

H = — ^tsC i+Sa C ia — t 2_^ C i+ri' o C in°\2 ~ * ' ^ >i + r i) 

+7^ 2^(Si " Sj+r? — -nin i+v ). (1) 

ir\a 

The first term contains hopping to first, second and third nearest neighbors (NN) with 
parameters t±, ti, t 3 respectively. The first NN of site i are labeled as i + rj. Eq. ([[]) 
is obtained from a standard canonical transformation of a Hubbard model with hoppings 
t\, £2, ts, if (complicated) terms smaller than t" = t 2 /U are neglected [16|. The difference 
between generalized t—J and strong-coupling Hubbard models is the meaning of the operator 
Cjcr, as explained below. The Hamiltonian can be written in terms of spinless fermions and 



spin- wave operators |B-ffl,^B|. We adopt the procedure and notation used by Martinez and 
Horsch [Q, slightly generalized to include second and third NN hoppings and the three- 
site term [|l(| : The sublattice A is defined as that of positive magnetization. The spins of 
sublattice B are rotated 180° around the x axis. In this way the Neel state is converted into a 
fully polarized ferromagnetic state, restoring the translational symmetry of the nonmagnetic 
state at the price of losing the conservation of spin. Then, the % operator is defined as a 
spinless holon creation operator h\, while c^ becomes a composite operator involving a local 
spin deviation aj. The result of both operations is the following representation: 

Cjj = h], Cii = h\ai , if j 6 A 

c iT = h\ai, Cii = h\ , if i € B. (2) 

In the exchange part (last term of Eq. (|TJ)) the fermion occupation numbers are averaged 
and the bosonic quadratic part is diagonalized by a standard canonical transformation: 



a q = u q a q - v q alg, (3) 

where u\ = v\ + 1 = 1/2 + l/(2i/,) , i/ g = (1 - 7j) 1 /2 j Uq > o, sgn(v g ) = sgn(7 g ), and ^ q = 
(cos q x + cos q y )/2. Retaining only linear terms in spin deviations for the rest of Eq. ([[]), 
the Hamiltonian becomes: 

H = E°j + Y^ u q a\a q + £ e k h[h k 
q k 

At 1 i 

+^L £ M(fc, ?)(/£/*_,«, + H.c), (4) 

where Sj is a constant, w g = 2 J u q , e k = (t 2 + 2(1 — x)t")e 2 (k) + (t 3 + (1 — x)t")ex(2k) and 
M(k,q) = (ug'jk-q + Vqlk), with ei(fc) = 47 fc and e 2 (k) = A cos /c x cos k y . In the present 
case, the doping x — 0. The constraint that at the same site there cannot be both a hole 
and a spin deviation is neglected since it does not affect the results for motion of a hole 
in a quantum antiferromagnet 0. The holon Green function Gh{k,uj) is obtained from the 
self-consistent solution of the following two equations: 

At, 
^K u) = -^ £ M 2 (k, q)G h (k -q,u- u q ) 

q 

G^ 1 (k, uj) = uj — e k — £(&, uj) + ie. (5) 

We have solved Eqs. (gj) in clusters of 16 x 16 and 20 x 20 sites. In order to obtain 
accurate values of the holon quasiparticle weight Zh, we have discretized the frequencies 
in intervals of Ao; = 10 _4 ti and have taken the small imaginary part e = 5Au. As an 
alternative method to that used by Liu and Manousakis ||], we have fitted the part of the 
spectral weight nearest to the quasiparticle peak by a sum of several Lorentzian functions. 
The resulting width of the quasiparticle peak was practically identical to 2e and from its 
integrated weight we determined Z^. We have verified that using this method there are 
practically no finite-size effects in our clusters. 

In the sudden approximation, the angle-resolved photoemission spectrum is proportional 
to the spectral density of states for Cu and O at wave vector k. These in turn are related 
to the imaginary part of the Green function for the generalized t — J operator c k(T or the 



generalized Hubbard operator Cka through a low-energy reduction procedure [12, Tq]. In linear 



order in 1/U, the well known procedure of the canonical transformation ||T4| , |T7|| applied to 
the generalized Hubbard model, in the subspace of no double occupancy, leads to: 

8 U 

Calling | 0) (| ipk)) the ground state of Eq. (f|) for the undoped (hole doped with wave 
vector k ) system, and using the Lehmann representation of the wave function, one realizes 
that while the holon quasiparticle weight is: 

Z h (k) =| (V* | h[ | 0) | 2 , (7) 

the weight for emitting a Hubbard electron is: 

Zg H (k) =| (V fc | ~c ka | 0) | 2 + | (ij k+Q | c ka | 0) | 2 , (8) 

where Q = (vr,7r), and | ipk) and | ipk+o) are the degenerate eigenstates of lowest energy of 
Eq. (f|) with a finite overlap with c ka | 0) . The corresponding result for the generalized 
t — J model ZgJ J (k) is obtained taking infinite U. Since Z c ^{k) = Z c ^{k) we restrict to spin 
up in the following. The states | ipk) can be constructed following the procedure used by 
Reiter ||. The only change in Eqs. 1 to 10 of Ref. ||, is that the quasiparticle energy Xk = 
Xk+Q is replaced by \k — e^ in Eqs. 3, 6 and 9, and by Xk — tk-q in Eq. 4 . Thus, writing 
explicitely only the terms with less than two spin-wave excitations we have: 

| Vfc) = A (k)h{ | 0) + ^- E 9 A,(k, q)h\_ q a\ | 0) + ..., (9) 

where: 

A t {k, q) = 4hM(k, q)G h {k -q,\ k - U q )A {k). (10) 

Using Eqs. (0) and (|B|) and retaining only terms lines in spin deviations we obtain: 

c iT = h\ - -i(l - x) J2 h\ +v at , if z G A 

U jj 

q t = h\di + — — [tt J2 h l+v + ^tshl+sidi+s ~ Oi)], i e B. (11) 



The most important correction of order 1/U is the first term between brackets in the second 
Eq. (PD and reflects the fact that in the ground state of the undoped Hubbard model, there 
is a finite double occupancy at sites B and an electron with spin up can be destroyed there, 
leaving a hole in one of its NN (this leads to the second term between brackets in Eqs. ( |P2|) 
and (JTJ)). 

Expressing Eqs. (|TT| ) m Fourier components, and using Z)ieA(B) e%kRi — (<5fc,o + 
e iQR '5 kyQ )N/2, we obtain: 



4 T = ~(1 + f(k))(hl + s A h\ +Q ) + ^= Y,(h{ +q - s A h{ +q+Q )[(l + g(k, q))a q - /(fc)o{], 

(12) 



where the phase sa = e l ® Ri with i e A, and 
f(k) = |(1 - x)e 1 (k), g(k, q) = l -^[t 2 {e 2 {k) - e 2 (k + q)) + t 3 (e 1 (2k) - e,(2k + 2q))]. 

(13) 

Using Eqs. ©, ©, ©, ®, © and (0) we obtain the desired result: 

-f^ = ^\l + f(k) + 1 ^j:M(k,q)G h (k~q 1 X k ~u Jq )[v q (l+g(k,q))-u q f(k)]\ 2 . 

(14) 

The sum is restricted to the magnetic Brillouin zone and the term with q = is excluded 
(there are no magnons with q = or q = Q in the | ipk) )• The weight Z^ J for the 
generalized t — J model operator Cj CT is given by Eq. ([14]) with the Hubbard perturbative 
corrections f(k) (first NN) and g(k, q) (second and third NN) set to zero. 

In Fig. 1 we compare the weight for the t — J model obtained by exact diagonalization 
Z^ D (k) in a square lattice of 20 sites [0 with our results Z^(k) for the 20 x 20 cluster at 



equivalent wave vectors. The comparison between exact results for square clusters of 16, 18, 
20 and 26 sites suggest that while the Z^ D (k) are nearly 20% larger for the 4x4 cluster, the 
finite size effects are of the order of 5% for larger clusters [I3|. The agreement between the 



exact Z l ^ D {k) and SCBA Z^(k) results is quite satisfactory. Note that the very small value 

6 



of Z tJ (Q) is a severe test to Eq. (|T4|), since it requires a near cancellation of the different 
terms. Instead, the "bare" SCBA result satisfies Z h {k) = Z h {k + Q) and cannot reproduce 
the shape of the exact results. 

With the confidence gained by the above comparison, we have calculated the generalized 
t — J and Hubbard weights for parameters which fit the observed quasiparticle dispersion 
Afc in Sr 2 Cu02Cl2 0. There are several choices of t 2 ,t 3 and t", including different signs 
of t" which produce nearly identical results. We took the parameters of Ref. ||. The 
resulting dispersion and weights are represented in Fig. 2. Compared with the parameters 
of Fig. 1, the effects of £2, £3 and t" are dramatic. They push the A& towards the incoherent 
part of the spectrum and reduce considerably the weights for the lowest A& (in the electron 
representation of Fig. 2). As a consequence, we could not detect quasiparticles near k = 0, 
Q or (it, 0) (Zh < 10~ 4 for these k). Therefore, the corresponding A& are not represented 
in Fig. 2. The weights for the generalized t — J and Hubbard models have significant 
differences: in contrast to the results for t 2 = £3 = t" = (not shown), Z^ J {k) is larger for 
k = (it/2 + e, it/ 2 + e) than for k = (n/2 — e, n/2 — e) with small e. Instead, Z^ a H (k), in 
agreement with experiment, is larger inside the non interacting Fermi surface. This effect is 
more noticeable for smaller values of U {t\/U = 0.1 was taken in Fig. 2) | |19| |. 

In summary, using the SCBA and related wave function, we have calculated the disper- 
sion and quasiparticle weight for removing a real electron in an undoped antiferromagnet 
described by a generalized t — J or a generalized Hubbard model in the strong coupling limit. 
The weight for the t — J model agrees very well with available exact results in sufficiently 
large clusters. While the generalized Hubbard can explain well both the measured disper- 
sion and weight of the quasiparticle in Sr 2 Cu02Cl2, the generalized t — J model, without 
mapping the electron operators, cannot. 

One of us (FL) is supported by the Consejo Nacional de Investigaciones Cientificas y 
Tecnicas (CONICET), Argentina. (AAA) is partially supported by CONICET. 

Note added: after submission of this manuscript we became aware of exact diagonalization 
results of the t — J model in a square cluster of 32 sites with periodic boundary conditions 



which has 9 non-equivalent wave vectors f20| . (See Fig. 3). The dispersion relation Afc agrees 
very well with the SCBA results except at the points k = (0,0), (7r/4, 7r/4) and (7r, 7r/2), 
where finite-size effects are obvious from the fact that A& ^ Xk+Q- Except at k = (0, 0) and 
(vr/4, 7r/4), where the position of A^ affects the quasiparticle weights, these weigths are in 
excellent agreement with our results using Eq. (|T4]). 
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FIGURE CAPTIONS 

Fig.l: Quasiparticle weight of the t — J model Z*£(k) calculated with the SCBA in a 20 x 20 
lattice for several wave vectors (triangles) , compared with exact diagonalization results in 
a square cluster of 20 sites Z^ D (k) |13|JT8|| (squares), and the spinless holon weight Zh(k) of 
the SCBA (circles). Parameters are t\ — 1, J = 0.3, t 2 = t 3 = t" = 0. 
Fig. 2: Top: quasiparticle dispersion in clusters of 16 x 16 (solid symbols) and 20 x 20 
sites (open symbols). Bottom: corresponding generalized t — J (squares) and generalized 
Hubbard (circles) quasiparticle weights. Parameters are: t\ = 0.35, t 2 = —0.12, t 3 = 

10 



0.08, J = 0.15, t" = J/4 and U = 3.5 

Fig. 3: Top: quasiparticle dispersion in the cluster of 16 x 16 (open squares) compared with 
the exact diagonalization results in a square cluster of 32 sites |2(J (solid triangles). Bottom: 
corresponding quasiparticle weights. 
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